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NATO  RTO  Working  Group  10  was  established  in  December  1998  to  examine  the 
technologies  of  plug  nozzles,  scramjet  propulsion  and  Computational  Fluid  Dynamics 
(CFD)  for  design  of  propelled  hypersonic  vehicles.  The  CFD  subgroup  evaluated  six 
topics:  boundary  layer  transition  and  instability,  real  gas  flows,  laminar  viscous-inviscid 
interactions,  shock-shock  interactions,  shock  wave  turbulent  boundary  layer  interactions 
(SWTBLI)  and  base  flows  with  and  without  plume  interaction.  This  paper  presents  a 
summary  of  the  evaluation  of  shock  wave  turbulent  boundary  layer  interactions.  Five 
configurations  were  considered:  2-D  compression  corner,  2-D  expansion-compression  cor¬ 
ner,  2-D  shock  impingement,  3-D  single  fin  and  3-D  double  fin.  Recent  Direct  Numerical 
Simulations  (DNS),  Large  Eddy  Simulations  (LES)  and  Reynolds-averaged  Navier-Stokes 
(RANS)  simulations  are  compared  with  experiment.  The  capabilities  and  limitations  are 
described,  and  future  research  needs  identified. 
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cf 

Skin  friction  coefficient 

P 

Ch 

Heat  transfer  coefficient 

aw 

Cp 

Pressure  coefficient 

c 

k 

Turbulence  kinetic  energy 

DNS 

L 

Length 

E 

i 

Turbulence  length  scale 

LES 

M 

Mach  number 

LES-D 

P 

Pressure 

LES-M 

Re 

Reynolds  number 

LES-X 

T 

Static  temperature 

LES-S 

U ,  V ,  w 

Cartesian  velocity  components 

m 

x:y:z 

Cartesian  coordinates 

RANS 

uT 

Wall  friction  velocity 

SA 

a 

Corner  angle 

sep 

S 

Boundary  layer  thickness 

ss 

K 

Von  Karman’s  constant 

SWTBLI 

P 

Density 

TA 
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Shear  stress 

Dynamic  molecular  viscosity  (=  up) 
Adiabatic  wall 
Characteristic  length 
Direct  Numerical  Simulation 
Experiment 

Large  Eddy  Simulation 
Dynamic  Smagorinsky  model  LES 
Monotone  Integrated  Large  Eddy  Simulation 
Mixed  scale  model  LES 
Static  Smagorinsky  model  LES 
Mean  value 

Reynolds-averaged  Navier-Stokes 

Spatial  accuracy 

Separation  length 

Converged  to  steady  state 

Shock  wave  turbulent  boundary  layer 

interaction 

Temporal  accuracy 

Throat  middle  line 

Van-Driest  transformation 

Wall 

Upstream 

Favre-fluctuating  component 
Normalized  by  vw/uT 
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Introduction 


Effective  design  of  supersonic  and  hypersonic  air 
vehicles  requires  accurate  simulation  methods  for  pre¬ 
dicting  aerothermodynamic  loads  (z.e.,  mean  and  fluc¬ 
tuating  surface  pressure,  skin  friction  and  heat  trans¬ 
fer).  Shock  wave  turbulent  boundary  layer  interaction 
(swtbli)  is  common  in  high  speed  flight,  and  can  sig¬ 
nificantly  affect  the  aerothermodynamic  loads.  The  re¬ 
cent  AGARD  Working  Group  18  (WG  18)  study  edited 
by  Knight  and  Degrez1  examined  the  capability  of  gen¬ 
eral1  Reynolds-averaged  Navier-Stokes  (rans  )  meth¬ 
ods  for  prediction  of  2-D  and  3-D  shock  wave  turbulent 
boundary  layer  interactions.  The  study  comprised  13 
configurations  and  37  computations  contributed  by  18 
participants.  The  recent  review  by  Zheltovodov3  also 
examined  the  capability  of  RANS  models  for  prediction 
of  swtbli.  Both  studies  concluded  that  general  RANS 
methods  accurately  predicted  mean  surface  pressure 
and  heat  transfer,  and  primary  separation  location  for 
weak  to  moderate11  2-D  and  3-D  shock  wave  turbulent 
boundary  layer  interactions.  For  strong  interactions, 
computed  surface  pressure  and  heat  transfer  distribu¬ 
tions,  and  location  of  primary  separation,  show  signif¬ 
icant  disagreement  with  experiment.  In  addition,  the 
RANS  calculations  failed  to  capture  the  high  level  of 
unsteadiness  in  the  shock  system  observed  in  the  ex¬ 
periment  and  the  appearance  of  secondary  separation, 
and  no  RANS  calculation  provided  a  prediction  of  rms 
fluctuating  surface  pressure  and  heat  transfer. 

RTO  Working  Group  10  (WG  10)  was  established  in 
December  1998  to  examine  three  specific  technologies 
for  propelled  hypersonic  flight:  plug  nozzles  (Sub¬ 
group  1),  scramjet  propulsion  (Subgroup  2)  and  Com¬ 
putational  Fluid  Dynamics  (Subgroup  3).  Six  specific 
topics  were  selected  for  evaluation  by  Subgroup  3  (Ta¬ 
ble  1)  and  specific  experiments  identified  for  CFD 
validation.  A  description  of  test  cases  was  published.4 
The  test  cases  are  also  described  in  Knight.5 

This  paper  presents  the  summary  report  of  Topic 
No.  5,  Shock  Wave  Turbulent  Boundary  Layer  Inter¬ 
actions.  Five  separate  configurations  are  considered. 
The  2-D  compression  corner  (Fig.  1)  and  expansion- 
compression  corner  (Fig.  2)  represent,  for  example, 
deflected  aerodynamic  control  surfaces.  The  2-D  shock 
impingement  (Fig.  3)  typifies  a  shock  boundary  layer 
interaction  within  an  inlet,  for  example.  The  3-D  sin¬ 
gle  fin  (Fig.  4)  is  typical  of  a  fin-fuselage  juncture.  The 
3-D  double  fin  (Fig.  5)  represents  a  sidewall  compres¬ 
sion  inlet.  All  five  configurations  have  been  extensively 
studied  and  a  substantial  set  of  experimental  data  is 
available  for  each  configuration. 

*By  general  rans  method  we  mean  a  rans  model  which 
has  not  been  modified  to  take  into  account  specific  details  of 
the  flowfield  structure.  Significant  progress  has  been  made  in 
development  of  engineering  rans  methods  which  incorporate 


Fig.  1  2-D  compression  corner 


A 

Fig.  2  2-D  expansion-compression  corner 


separation  bubble 

Fig.  3  2-D  shock  impingement 


Fig.  4  3-D  single  fin 


Fig.  5  3-D  double  fin 


Contributions  of  recent111  numerical  simulations  for 

configuration-dependent  flow  physics.  An  example  is  Panaras.2 

11  In  this  context,  “weak”  and  “moderate”  refer  to  the  overall 
static  pressure  rise. 

“Simulations  were  selected  from  1997  onward,  since  the  re¬ 
view  of  Knight  and  Degrez1  for  AGARD  WG  18  was  completed 
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Table  1  Topics  and  Team  Leaders 


No, 

1 

2 

3 

4 

5 

6 


Topic 

Boundary  Layer  Instability  and  Transition 
Real  Gas  Flows 

Laminar  Viscous-Inviscid  Interactions 
Shock-Shock  Interactions 

Shock  Wave-Turbulent  Boundary  Layer  Interactions 
Base  Flows  with  and  without  Plume  Interactions 


Team  Leader (s) 

Steven  Schneider  (Purdue  Universty) 

Graham  Candler  (University  of  Minnesota) 
Michael  Holden  (CUBRC) 

Steven  Walker  and  John  Schmisseur  (AFOSR) 
Doyle  Knight  (Rutgers  University) 

Peter  Bakker (Technical  University  of  Delft) 
Phillipe  Reijasse  (ONERA) 


the  five  cases  were  solicited  through  the  RTO  WG  10 
membership  and  the  RTO  WG  10  website.6  The  stud¬ 
ies  include  Direct  Numerical  Simulation  (dns),  Large 
Eddy  Simulation  (les)  and  Reynolds-averaged  Navier 
Stokes  (rans).  The  computed  results  are  compared 
herein  with  the  experimental  data,  and  the  overall  ca¬ 
pabilities  of  the  CFD  methods  are  assessed.  Emphasis 
is  placed  on  validation  of  the  numerical  methods,  and 
therefore,  in  the  interest  of  space,  the  description  of 
the  flowfield  structures  is  intentionally  limited. 

2-D  Compression  Corner 

The  2-D  supersonic  compression  corner  (Fig.  6)  is  a 
standard  CFD  validation  case  for  shock  wave  turbulent 
boundary  layer  interactions.  An  equilibrium  super¬ 
sonic  turbulent  boundary  layer  approaches  a  compres¬ 
sion  corner  defined  by  angle  a.  The  deflection  of  the 
flow  by  the  corner  generates  a  shock  system.  For 
sufficiently  large  pressure  rise,  the  boundary  layer  sep¬ 
arates  and  a  A-shock  forms.  The  shock  system  displays 
significant  unsteadiness.7-10 


Fig.  6  2-D  compression  corner 


where  S  is  the  upstream  boundary  layer  thickness, 
is  the  upstream  freest  ream  Mach  number,  p2  is  the 
downstream  inviscid  pressure,  and  ppi  is  the  plateau 
pressure  computed  according  to  the  empirical  formula 
by  Zukoski15 


Ppi  — Poo  {y2MOQ  H- 1)  (2) 

The  solid  lines  define  the  limits  of  a  large  experimen¬ 
tal  database  for  supersonic  compression  corners  over 
a  wide  range  of  Mach  numbers.  The  rans  calcula¬ 
tions  are  listed  in  Table  2.  It  is  evident  that  rans 
computations  in  general  do  not  accurately  predict  the 
separation  length. 


An  extensive  number  of  analytical,  computational, 
and  experimental  studies  of  the  supersonic  compres¬ 
sion  corner  have  been  conducted.  Recent  reviews  in¬ 
clude  Zheltovodov3  and  Dolling.9,10  Earlier  reviews 
include  Green11  and  Delery.12 

General  rans  methods  have  failed  to  accurately 
predict  the  flowfield  structure  of  the  supersonic  com¬ 
pression  corner.  Fig.  7  displays  the  computed  mean 
separation  lengthlv,  normalized  by  the  characteristic 
length  Lc  proposed  by  Zheltovodov  and  Schiilein13,14 


r  /  \  3.1 

_  0  I  P2  \ 

C_Z  \Ppi ) 


(1) 


in  early  1997. 

lvThe  mean  separation  length  is  the  straight  line  distance 
between  the  mean  separation  and  attachment  points. 


Fig.  7  Separation  length  for  rans 

In  the  following  sections,  DNS  and  les  results  are 
presented  for  several  configurations.  Mean  quantities 
represent  combined  time  and  spanwise  averages. 

Case  1:  Mach  3,  a  =  8° 

Urbin  et  al 22  simulated  a  Mach  3,  a  =  8°  compres¬ 
sion  corner  at  Res  =  2  x  104  and  adiabatic  wall  con¬ 
ditions.  Experimental  data  of  Zheltovodov  et  a/23-26 
is  available  for  a  =  8°  at  Reynolds  numbers  Res  = 
7.5  x  104  and  1.1  x  105.  The  flow  conditions  are  summa¬ 
rized  in  Table  3.  The  inviscid  pressure  ratio  across  the 
shock  P2/P1  =  1.8.  As  discussed  in  Zheltovodov,3  this 
configuration  corresponds  to  an  incipient  intermittent 
separation.  Although  the  mean  flowfield  is  not  sepa¬ 
rated,  instantaneous  reversed  flow  can  exist  for  a  short 


Table  2  RANS  calculations  of  2-D  Compressionn  Corner 


Number 

Ref. 

Model 

a 

Mach 

Res  x  10“6 

Rsep/ $ 

Rsep/ Tc 

1 

Shang16 

Equilibrium 

25° 

2.96 

0.1 

5.0 

15.3 

2 

Shang16 

Frozen 

25° 

2.96 

0.1 

13.0 

39.7 

3 

Shang16 

Relaxation 

25° 

2.96 

0.1 

10 

30.5 

4 

Horstman17 

Baseline 

to 

o 

o 

2.85 

1.65 

1.45 

8.7 

5 

Horstman17 

Relaxation 

to 

o 

o 

2.85 

1.65 

2.69 

16.1 

6 

Horstman17 

Pressure  Gradient 

to 

o 

o 

2.85 

1.65 

0.81 

4.87 

7 

Horstman17 

Relaxation + 

to 

o 

o 

2.85 

1.65 

1.62 

9.75 

Pressure  Gradient 

8 

Horstman17 

Baseline 

to 

o 

2.85 

1.33 

4.0 

14.5 

9 

Horstman17 

Relaxation 

to 

o 

2.85 

1.33 

4.7 

17.1 

10 

Horstman17 

Pressure  Gradient 

to 

o 

2.85 

1.33 

2.74 

9.96 

11 

Horstman17 

Relaxation+ 

to 

o 

2.85 

1.33 

3.6 

13.1 

Pressure  Gradient 

12 

Visbal, Knight18 

B-L 

16° 

2.9 

1.6 

1.0 

14.4 

13 

Visbal,  Knight18 

Modified  B-L 

16° 

2.9 

1.6 

0.2 

2.89 

14 

Visbal, Knight18 

Relaxation 

16° 

2.9 

1.6 

0.6 

8.66 

15 

Visbal, Knight18 

Constant  Ymax 

16° 

2.9 

1.6 

0.2 

2.89 

16 

Visbal, Knight18 

B-L 

to 

o 

o 

2.9 

1.6 

1.6 

11.5 

17 

Visbal, Knight18 

Modified  B-L 

to 

o 

o 

2.9 

1.6 

0.9 

6.46 

18 

Visbal, Knight18 

Relaxation 

to 

o 

o 

2.9 

1.6 

1.8 

12.9 

19 

Visbal, Knight18 

Modified  B-L 

to 

o 

2.9 

1.6 

2.05 

7.25 

20 

Visbal, Knight18 

Relaxation 

to 

o 

2.9 

1.6 

2.95 

10.4 

21 

Ong, Knight19 

B-L 

16° 

1.96 

0.25 

2.15 

10.1 

22 

Ong, Knight19 

B-L 

16° 

2.83 

1.6 

0.6 

2.81 

23 

Ong, Knight19 

B-L 

to 

o 

o 

2.83 

1.6 

1.75 

8.19 

24 

Horstman20 

k  —  e 

25° 

3.0 

0.133 

4.89 

16.3 

25 

Borisov  et  al 21 

k  —  oo 

25° 

3.0 

0.139 

7.38 

24.8 

time  in  the  vicinity  of  the  shock  wave  and  near  the  wall 
due  to  the  unsteadiness  of  the  interaction  process. 


Table  3 

Reference 
Urbin  et  al 22 
Zheltovodov  et  al23~ 25 


Flow  Conditions 

Data  Moo  Re§  x  10“ 4 

les  3.0  2.0 

E  2.9  7.5-11 


Details  of  the  computation  are  presented  in  Table  4 
where  SA  and  TA  indicate  the  spatial  and  temporal 
accuracy.  An  unstructured  grid  LES  algorithm  with 
a  static  Smagorinsky  model  was  utilized.  The  inflow 
conditions  were  obtained  from  a  separate  computation 
of  a  flat  plate  adiabatic  boundary  layer  at  Mach  3.  The 
inflow  velocity  profile  satisfied  the  compressible  Law  of 
the  Wall.26  The  inflow  skin  friction  coefficient  agreed 
with  the  empirical  value  obtained  from  the  compress¬ 
ible  Law  of  the  Wall  and  Wake  to  within  6%,  and  the 
inflow  adiabatic  wall  temperature  was  within  3%  of  the 
theoretical  value  obtained  from  the  empirical  formula 


1  + 


(7-1) 


PrtmMi 


(3) 


where  Prt m  =  0.89  is  the  mean  turbulent  Prandtl 
number. 


Table  4  Details  of  Computations 

Reference  Type  Cells  SA  ta 

Urbin  et  al 22  les  1.6  x  106  2nd  2nd 

An  instantaneous  image  of  the  shock  wave  (as  visu¬ 
alized  by  the  isosurface  p/poo  =  1.25)  and  the  stream- 
wise  velocity  at  y+  =  10  are  shown  in  Fig.  8.  Regions 
of  negative  velocity  (instantaneous  separation)  are  vis¬ 
ible  in  qualitative  agreement  with  the  observations  of 
Zheltovodov.3 

The  mean  wall  pressure  distribution  is  shown  in 
Fig.  9.  The  experimental  profiles  for  the  two  different 
Reynolds  numbers  are  essentially  identical.  The  com¬ 
puted  profile  shows  close  agreement  with  experiment. 
The  mean  velocity  profiles  upstream  and  downstream 
of  the  corner  are  displayed  in  Fig.  10.  The  experi¬ 
mental  data  is  for  Res  =  1.1  x  105.  The  upstream 
profile  is  3.3<5  ahead  of  the  corner  and  in  the  undis¬ 
turbed  flat  plate  boundary  layer.  The  downstream 
profile  is  3.5 5  from  the  corner  and  located  close  to  the 
position  where  the  wall  pressure  reaches  its  asymp¬ 
totic  value.  The  computed  results  are  in  close  agree¬ 
ment  with  experiment.  The  mean  static  temperature 
profiles  at  the  same  locations  are  shown  in  Fig.  11. 
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Upstream  of  the  interaction  the  computed  and  exper¬ 
imental  profiles  are  in  close  agreement.  Downstream 
of  the  interaction,  the  profiles  agree  except  within  the 
region  0.05  <y/S<  0.25  where  the  experimental  value 
exceeds  the  prediction. 


Fig.  8  Instantaneous  velocity  and  shock 


Fig.  9  Mean  wall  pressure 


Fig.  10  Mean  velocity 


The  Reynolds  streamwise  and  shear  stresses  are  dis¬ 
played  in  Figs.  12  and  13,  respectively,  where  the 
stresses  have  been  normalized  by  the  upstream  wall 
shear  stress  rw.  The  experimental  data  of  Muck  et 
al 27  is  shown.  The  experimental  Reynolds  number 
Res  =  1.6  x  106  is  a  factor  of  eighty  greater  than  the 
les.  However,  scaling  of  the  Reynolds  shear  stress  in 
a  flat  plate  boundary  layer  by  the  wall  shear  stress 
rw  has  been  shown  to  collapse  the  experimental  data 
for  a  range  of  Mach  and  Reynolds  numbers.28  Thus, 


Fig.  11  Mean  static  temperature 

the  comparison  of  the  les  and  experiment  is  justi¬ 
fied  if  only  for  qualitative  comparison.  The  effect  of 
the  shock  boundary  layer  interaction  is  to  increase 
the  level  of  Reynolds  normal  pu"u 11  and  shear  pu"v 11 
stresses.  The  predicted  values  are  in  good  agreement 
with  experiment  except  for  y/S  <  0.1. 

In  summary,  the  following  conclusions  can  be  made: 

•  les  of  the  Mach  3,  a  =  8°  compression  cor¬ 
ner  shows  good  agreement  with  the  experimental 
mean  surface  pressure  distribution  and  qualita¬ 
tively  correct  modification  of  the  mean  velocity 
and  turbulence  profiles  by  the  shock  wave 

•  The  mean  flowfield  is  unseparated  (and  conse¬ 
quently,  this  case  does  not  represent  a  particularly 
difficult  case  for  les),  but  the  instantaneous  flow- 
field  indicates  an  incipient  intermittent  separation 
regime  in  accordance  with  the  generalization  of 
experimental  data  and  analysis  by  Zheltovodov.3 


Fig.  12  Reynolds  streamwise  stress 


Case  2:  Mach  3,  a  =  18° 

Adams29-31  and  Rizzetta  et  al32  simulated  a 
Mach  3,  a  =  18°  compression  corner  at  Res  =  2.1  x  104 
and  adiabatic  wall  conditions.  No  experimental  data  is 
available  for  this  configuration.  The  inviscid  pressure 
ratio  across  the  shock  is  p2/pi  =  3.37. 

Details  of  the  computations  are  presented  in  Table  6. 
Adams  performed  a  DNS.  Rizzetta  et  al  used  three 
different  approaches,  namely,  a  DNS,  an  les  using  a 


Fig.  13  Reynolds  shear  stress 


Table  5  Flow  Conditions 

Reference  Data  Res  x  10“ 4 

Adams29-31  dns  3  2.1 

Rizzetta  et  al 32  dns  ,  les  3  2.1 

static  Smagorinsky  model  for  the  SGS  stresses  (les- 
s),  and  an  les  using  a  dynamic  (Germano  et  al33) 
Smagorinksy  model  (les-d).  The  inflow  conditions 
were  obtained  in  each  instance  from  a  separate  com¬ 
putation  of  a  flat  plate  adiabatic  boundary  layer  at 
Mach  3  (Fig.  14). 

Table  6  Details  of  Computations 


Reference 

Type 

Grid 

SA 

TA 

Rizzetta  et  al 32 

DNS 

421x151x81 

6th 

2nd 

LES-S 

421x151x81 

6th 

2nd 

les-d 

421x151x81 

6th 

2nd 

Adams29-31 

DNS 

1000x180x80 

5th 

3rd 

flat-plate  domain 


l  l 

L 

half-channel  compression-ramp 

domain  domain 

Fig.  14  Generation  of  inflow  conditions 

Although  the  Mach  number  and  Reynolds  number 
were  nominally  identical,  the  inflow  profiles  differ  sig¬ 
nificantly.  Figs.  15  and  16  display  the  mean  and  rms 
fluctuating  streamwise  velocity  at  the  inflow  boundary 
(“present”  indicates  the  results  of  Rizzetta  et  al).  Ma¬ 
jor  differences  are  evident,  and  have  a  significant  effect 
on  the  computed  shock  boundary  layer  interaction. 

The  computed  mean  surface  pressure  coefficient 
cp  —  ip— Poo)l\PooU%o  is  shown  in  Fig.  17.  The  size 
of  the  interaction  region,  defined  by  the  location  of 
the  initial  rise  in  cp  and  intersection  of  cv  with  the 
downstream  inviscid  value  (cp  =  0.376),  is  substan¬ 
tially  larger  in  the  simulation  of  Adams  compared  to 
Rizzetta.  This  is  attributable  to  the  differences  in  the 
inflow  velocity  profile  (Fig.  15),  as  the  momentum  flux 
pu  is  lower  near  the  wall  in  Adams’  profile  (see  also 


Fig.  15  Inflow  mean  velocity 


Fig.  16  Inflow  rms  streamwise  velocity 

Fig.  20).  The  dns  and  les  profiles  of  Rizzetta  are  vir¬ 
tually  identical,  implying  the  SGS  model  has  negligible 
effect. 


Fig.  17  Mean  surface  pressure 

The  computed  mean  skin  friction  coefficient  Cf  = 
rw / \ Poo is  shown  in  Fig.  18.  The  upstream  skin 
friction  level  is  virtually  identical  for  both  Adams  and 
Rizzetta  (see  Fig.  15  which  indicates  the  velocity  pro¬ 
files  are  essentially  the  same  within  the  viscous  sub¬ 
layer).  The  separation  length  predicted  by  Adams  is 
approximately  130%  larger  than  Rizzetta.  This  is  at¬ 
tributable  to  the  lower  momentum  flux  in  the  incoming 
boundary  layer  (Fig.  15).  The  dns  and  les  profiles  of 
Rizzetta  are  virtually  identical,  again  indicating  the 
SGS  model  has  negligible  effect. 

The  mean  velocity  and  static  temperature  are  dis¬ 
played  in  Figs.  19  and  20,  respectively.  Station  Nos.  1, 
2  and  3  are  located  upstream  of  the  interaction  re- 


Fig.  18  Mean  skin  friction 

gion,  at  the  corner,  and  downstream  of  the  corner, 
respectively.  The  mean  velocity  profiles  are  displayed 
using  the  Van  Driest34  transformation.  All  profiles  ex¬ 
hibit  agreement  with  the  Law  of  the  Wall  logarithmic 
profile  in  the  region  20  <  y+  <  100  (Station  No.  1). 
Downstream  of  the  corner,  the  simulations  by  Rizzetta 
show  recovery  of  the  near  wall  region  to  the  equilib¬ 
rium  Law  of  the  Wall,  while  the  simulation  by  Adams 
does  not.  Significant  differences  in  the  mean  temper¬ 
ature  between  Adams  and  Rizzetta  are  evident  both 
upstream,  at  and  downstream  of  the  corner.  These  are 
associated  with  differences  in  the  mean  velocity,  since 
the  total  temperature  is  approximately  constant  across 
the  boundary  layer.  For  both  the  velocity  and  temper¬ 
ature  profiles,  the  dns  and  les  results  of  Rizzetta  are 
virtually  identical. 


Fig.  20  Mean  static  temperature 


Spanwise  turbulent  kinetic  energy  spectra  at  two 
separate  positions  within  the  boundary  layer  and  three 
streamwise  locations  are  shown  in  Figs.  21  and  22 
where  the  location  is  defined  in  terms  of  the  inflow 
boundary  layer  condition.  All  curves  have  been  nor¬ 
malized  to  unity  at  &3  =  1.  The  results  of  Rizzetta 
are  consistently  higher  than  Adams,  implying  a  higher 
level  of  turbulent  kinetic  energy.  This  result  is  con¬ 
firmed  by  comparison  of  profiles  for  the  streamwise, 
normal  vertical  and  spanwise  rms  fluctuating  veloci¬ 
ties.32  The  profiles  of  Rizzetta  display  the  emergence 
of  an  inertial  subrange  (E(k)  oc  &-5/3).  The  inertial 
subrange  properly  occurs  at  wavenumbers  above  the 
peak  in  the  energy  spectrum.35  Adams’  results  do  not 
display  an  inertial  subrange.  The  spectra  of  Rizzetta 
at  high  wavenumber  show  higher  energy  for  the  LES 
than  dns  calculations.  This  is  an  unexplained  result, 
since  it  would  be  expected  that  the  SGS  model  would 
dissipate  more  energy  at  higher  wavenumber  than  the 
DNS. 


- DNS 


Fig.  21  Spanwise  TKE  spectra  at  y+  =  7.3 (y/S  = 
0.034) 


- DNS 


Fig.  22  Spanwise  TKE  spectra  at  y+  —  132 (y/S  — 
0.61) 

The  computed  separation  length,  normalized  using 
(1),  is  shown  in  Fig.  23.  The  computed  separation 
lengths  of  Rizzetta  are  below  the  linear  extrapolation 
of  the  lower  bound  of  the  experimental  band,  while  the 
computed  separation  length  of  Adams  lies  above  the 
upper  bound. 

In  summary,  the  following  conclusions  may  be  made: 
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Table  8 


Details  of  Computations 


Fig.  23  Separation  length  for  DNS  and  LES 

•  DNS  of  shock  wave  turbulent  boundary  layer  in¬ 
teractions  is  sensitive  to  the  inflow  profiles 

•  The  DNS  and  LES  of  Rizzetta  are  virtually  iden¬ 
tical  for  all  profiles.  This  implies  that  the  SGS 
model  has  a  negligible  effect. 

•  The  computed  mean  separation  lengths  are  out¬ 
side  the  linear  extrapolation  of  the  experimental 
data 

Case  3:  Mach  3,  a  =  24° 

Rizzetta  and  Visbal36  simulated  a  Mach  3,  a  =  24° 
compression  corner  at  Res  =  2.1  x  104  and  adia¬ 
batic  wall  conditionsv.  Experimental  data  is  available 
from  Settles  et  a/,37  Dolling  and  Murphy38  and  Smits 
and  Muck39  as  indicated  in  Table  7V1.  However,  the 
Reynolds  number  for  the  experiment  is  significantly 
higher  than  the  simulation  by  a  factor  of  40  to  77. 
Therefore,  the  comparison  between  simulation  and  ex¬ 
periment  is  qualitative.  The  inviscid  pressure  ratio 
across  the  shock  is  p2/pi  =  4.68. 


Reference  Type  Grid  SA  ta 

Rizzetta  et  al 36  les  421  x  151  x  181  6th  2nd 

(Fig.  24),  streamwise  Reynolds  stress  (Fig.  25)  and 
Reynolds  shear  stress  (Fig.  26).  The  computed  inflow 
mean  velocity  is  lower  than  the  experiment  near  the 
wall  as  expected. 


Fig.  24  Inflow  velocity 


Fig.  25  Inflow  Reynolds  streamwise  stress 


Table  7  Flow  Conditions 


Reference 
Rizzetta  et  al 36 
Settles  et  al37 
Dolling  et  al38 
Smits  et  al 39 


Data  Moo 

les  3.0 

E  2.84  -  2 

E  2.90  -  2 

E  2.79  -  2 


Res  x  1(T6 

.021 

.87 

1.5 -1.6 

.95 

0.81  -  1.4 

.87 

1.6 

Details  of  the  computation™  are  presented  in  Ta¬ 
ble  8.  A  structured  grid  LES  algorithm  using  a  dy¬ 
namic  Smagorinsky  model  was  utilized. 

The  effect  of  the  difference  in  Reynolds  number  is 
evident  in  the  inflow  profiles  of  streamwise  velocity 

vRizzetta  and  Visbal36  also  simulated  a  —  8°,  16°  and  20° 
at  Res  —  2.1  x  104. 

V1  Please  note  that  the  reference  numbers  in  Figs.  24  to  32 
correspond  to  the  references  in  Rizzetta  and  Visbal.36 

V11A  coarse  grid  simulation  (211  x  76  x  81)  was  also  performed 
which  demonstrated  that  the  fine  grid  solution  for  the  surface 
pressure  was  grid  converged. 


- LES 

°  experiment  9 


~Pm(Ut  Un)m 


0.06 


Fig.  26  Inflow  Reynolds  shear  stress 

The  mean  wall  pressure  is  shown  in  Fig.  27.  The 
upstream  influence,  plateau  and  recovery  regions  are 
accurately  predicted,  notwithstanding  the  difference 
in  Reynolds  number  between  the  simulation  and  ex- 


periment.  The  standard  deviation  ap  of  the  surface 
pressure  fluctuations,  normalized  by  the  local  mean 
wall  pressure  pw>  is  displayed  in  Fig.  28.  The  exper¬ 
imental38  peak  value  is  associated  with  the  fluctuat¬ 
ing  motion  of  the  separation  shock.  The  simulation 
significantly  underestimates  the  peak  value,  and  over¬ 
estimates  the  level  of  fluctuations  downstream  of  the 
corner.  Similarly,  Fig.  29  indicates  that  the  simulation 
underestimates  the  peak  value  of  skewness  S  =  p'3  /a3 
near  separation.  These  discrepancies  can  be  under¬ 
stood  by  considering  the  time  history  of  the  surface 
pressure  fluctuations  in  the  experiment  and  simula¬ 
tion.  Fig.  30  shows  a  time  history  of  the  surface 
pressure  at  midspan  and  x/S  =  —2.1  which  is  near 
the  location  of  the  peak  values  of  the  standard  devia¬ 
tion  (Fig.  28)  and  skewness  (Fig.  29).  The  experiment 
displays  a  high  intermittency  of  the  pressure  signal  as¬ 
sociated  with  low  frequency  motion  of  the  separation 
shock.  The  simulation  displays  no  low  frequency  com¬ 
ponent. 


Fig.  27  Mean  wall  pressure 


Fig.  28  Standard  deviation  of  wall  pressure 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  31. 
The  position  of  separation  and  reattachment  (as  mea¬ 
sured  from  the  corner  at  x  =  0)  are  significantly 
overestimated,  and  the  recovery  of  the  skin  friction 
downstream  of  attachment  substantially  underpredicts 
the  experiment. 

The  computed  separation  length,  normalized  using 
(1),  is  shown  in  Fig.  23.  The  computed  separation 


Fig.  29  Skewness  of  wall  pressure 


t/ 6, 


Fig.  30  Time  history  of  wall  pressure 


Fig.  31  Mean  skin  friction 

length  is  at  the  lower  limit  of  a  linear  extrapolation  of 
the  experimental  band. 

Velocity  profiles  at  two  locations  downstream  of  the 
corner  at  x/5  —  1.33  and  2.65  are  displayed  in  Fig.  32. 
The  underestimate  in  the  velocity  in  the  near  wall  re¬ 
gion  is  consistent  with  the  comparison  of  computed 
and  experimental  upstream  profiles  (Fig.  24). 

In  summary,  the  following  conclusions  can  be  made: 

•  The  mean  surface  pressure  shows  good  agree¬ 
ment  with  experiment.  This  may  be  fortuitous, 
however,  in  view  of  the  significant  difference  in 
Reynolds  number  between  the  simulation  and  ex¬ 
periment,  and  the  substantial  differences  in  sur¬ 
face  pressure  fluctuation  statistics  in  the  vicinity 
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Fig.  32  Velocity  profiles 

of  separation. 

•  The  computed  separation  length  is  at  the  lower 
limit  of  the  experimental  data  correlation. 

•  The  near  wall  velocity  profiles  downstream  of  the 
corner  are  significantly  under  predicted,  likely  due 
to  the  difference  in  inflow  profiles  associated  with 
the  disparity  in  Reynolds  numbers. 

Case  4:  Mach  3,  a  =  25° 

Yan  et  al 40  simulated  a  Mach  2.88,  a  =  25°  com¬ 
pression  corner  at  Res  =  2  x  104  and  adiabatic 
wall  conditions.  Experimental  data  of  Zheltovodov 
et  a/3,23-25,41  is  available  for  a  =  25°  at  Reynolds 
number  Res  =  6.4  x  104.  The  flow  conditions  are 
summarized  in  Table  9.  The  inviscid  pressure  ratio 
across  the  shock  p2/pi  =  4.7.  The  mean  flowfield  is 
separated  at  the  corner. 

Table  9  Flow  Conditions 

Reference  Data  Res  x  10“ 4 

Yan  et  al 40  les  2.88  2.0 

Zheltovodov  et  al3,23~25,41  E  2.88  6.4 

Details  of  the  computation  are  presented  in  Ta¬ 
ble  10.  An  unstructured  grid  LES  algorithm  using  the 
miles  model42  was  utilized.  In  the  miles  model,  the 
inherent  dissipation  in  the  numerical  algorithm  con¬ 
stitutes  the  SGS  model.  The  inflow  velocity  profile 
satisfied  the  Law  of  the  Wall.  The  inflow  skin  friction 
coefficient  agreed  with  the  empirical  value  obtained 
from  the  compressible  Law  of  the  Wall  and  Wake  to 
within  6%.  The  inflow  adiabatic  wall  temperature  was 
within  3%  of  the  theoretical  value  obtained  from  (3). 

Table  10  Details  of  Computations 

Reference  Type  Cells  SA  ta 

Yan  et  al 40  les  2.0  x  106  2nd  2nd 

The  mean  wall  pressure  is  shown  in  Fig.  33vm.  The 
upstream  influence  (i.e.:  the  location  of  the  initial 

vlllThe  reference  number  for  the  experiment  in  Figs.  33  and  34 
correspond  to  the  references  listed  in  Yan  et  al.40 


pressure  rise)  is  accurately  predicted;  however,  the 
computed  pressure  distribution  does  not  show  the  dis¬ 
tinct  plateau  evident  in  the  experiment,  and  recovers 
to  the  downstream  inviscid  pressure  more  rapidly  than 
the  experiment. 


Fig.  33  Mean  wall  pressure 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  34. 
The  initial  drop  in  skin  friction  associated  with  the 
adverse  pressure  gradient  is  reasonably  predicted.  The 
predicted  recovery  of  the  skin  friction  downstream  of 
reattachment  is  more  rapid  than  in  the  experiment, 
similar  to  the  results  for  the  surface  pressure. 


Fig.  34  Mean  skin  friction  coefficient 

The  computed  mean  separation  length  is  shown  in 
Fig.  23.  The  value  lies  within  a  linear  extrapolation  of 
the  band  of  experimental  data. 

Fig.  35  displays  the  instantaneous  shock  wave  struc¬ 
ture  defined  by  an  isosurfaces  of  static  pressure  at 
p/Poo  =  1-4  and  2.0,  together  with  the  contours  of  the 
streamwise  velocity  at  y+  =  7.  The  A  shock  is  consis¬ 
tent  with  the  experimental  Schlieren  image  (Fig.  36). 

In  summary,  the  following  conclusions  can  be  made: 

•  The  les  does  not  accurately  reproduce  the  mean 
surface  pressure.  This  may  be  attributable  to  the 
difference  in  Reynolds  numbers  between  the  les 
and  experiment. 


Fig.  36  Schlieren  image 

•  The  LES  predicts  a  mean  separation  length  con¬ 
sistent  with  a  linear  extrapolation  of  the  experi¬ 
mental  data. 

2-D  Expansion- Compression  Corner 

The  expansion-compression  corner  (Fig.  37)  pro¬ 
vides  a  test  of  modeling  of  shock  wave  turbulent 
boundary  layer  interaction  wherein  the  incoming 
boundary  layer  is  non-equilibrium.  An  initially  equi¬ 
librium  supersonic  turbulent  boundary  layer  is  ex¬ 
panded  through  an  angle  a  and  subsequently  com¬ 
pressed  through  the  same  angle.  The  expansion  re¬ 
duces  the  turbulence  kinetic  energy  and  increases  the 
mean  kinetic  energy  within  the  boundary  layer. 

-Mqq  / 


number  downstream  of  the  expansion  is  M2  =  4.59. 
Across  the  subsequent  compression,  the  pressure  ra¬ 
tio  is  P3/P2  =  8.60  and  the  downstream  Mach  number 
M3  =  2.43. 


Table  11  Flow  Conditions 

Reference  Data  Res  x  10“ 4 

Yan  et  al 43  LES  2.88  2.0 

Zheltovodov13,14,25,44-48  E  2.88  4.4  -  50. 


Details  of  the  computation  are  presented  in  Ta¬ 
ble  12.  An  unstructured  grid  LES  algorithm  using  the 
miles  model42  was  utilized.  The  inflow  profile  is  the 
same  as  Case  1  of  the  2-D  compression  corner. 

Table  12  Details  of  Computations 

Reference  Type  Cells  SA  TA 

Yan  et  al43  les  2.4  x  106  2nd  2nd 


The  computed  flowfield  structure  is  shown  in 
Figs.  38  and  Fig.  39  which  display  the  mean  static 
pressure  and  streamlines  (s  is  separation,  A  is  attach¬ 
ment).  The  flow  expands  around  the  first  corner, 
and  recompresses  at  the  second  corner.  The  shock 
wave  turbulent  boundary  layer  at  the  second  corner  is 
sufficiently  strong  to  separate  the  boundary  layer  as 
evident  in  Fig.  39.  The  flowfield  structure  is  in  good 
agreement  with  the  results  of  Zheltovodov  and  his  col¬ 
leagues3,  13, 14,25,45, 47, 48  which  are  shown  qualitatively 
in  Fig.  37. 
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Fig.  37  Expansion-compression  corner 

In  the  following  section,  an  LES  of  the  expansion- 
compression  corner  is  presented.  Mean  quantities  rep¬ 
resent  combined  time  and  spanwise  averages. 

Case  1:  Mach  3,  a  =  25° 

Yan  et  al 43  simulated  a  Mach  2.88  expansion  com¬ 
pression  corner  at  Res  =  2  x  104  and  adiabatic 
wall  conditions.  Experimental  data  of  Zheltovodov  et 
a^i3,i4,25,44-48  ayaiiable  for  a  =  25°  at  =  2.9 

and  Res  =  4.1  x  104  to  1.95  x  105.  The  flow  conditions 
are  shown  in  Table  11.  The  inviscid  pressure  ratio 
across  the  expansion  is  p2/pi  =  0.098,  and  the  Mach 


Fig.  38  Mean  static  pressure 


Fig.  39  Mean  streamlines 


Mean  velocity  profiles  in  the  ^-direction  are  shown 
in  Fig.  40  at  x  =  25  and  x  =  65,  where  x  is  mea¬ 
sured  from  the  inflow  along  the  direction  of  the  inflow 
freestream  velocity  (Fig.  38).  The  abscissa  is  the  com¬ 
ponent  of  velocity  locally  parallel  to  the  wall,  and 
the  ordinate  is  the  distance  measured  normal  to  the 
wall.  The  first  profile  is  upstream  of  the  expansion 
corner  (which  is  located  at  x  =  45).  The  second  pro¬ 
file  is  downstream  of  the  expansion  fan  and  upstream 
of  the  separation  point.  The  computed  mean  velocity 
profile  at  the  first  location  is  slightly  fuller  than  the  ex¬ 
periment.  This  is  consistent  with  the  experimentally 
observed  dependence  of  the  exponent  n  in  the  power- 
law  U/Uoq  =  ( y/5 )1/n  on  the  Reynolds  number.  The 
second  profile  shows  a  significant  acceleration  of  the 
flow  in  the  outer  portion  of  the  boundary  layer  due  to 
the  expansion. 


The  surface  pressure  is  displayed  in  Fig.  41.  The 
surface  pressure  drops  rapidly  at  the  first  corner.  The 
shock  boundary  layer  interaction  at  the  second  cor¬ 
ner  causes  separation  of  the  boundary  layer  and  a 
concomitant  plateau  in  the  surface  pressure.49  The 
experiments  display  an  increase  in  the  size  of  the  pres¬ 
sure  plateau  region  with  decreasing  Reynolds  number. 
The  computed  results  for  Res  =  2  x  104  are  in  good 
agreement  with  the  experimental  data  at  the  lowest 
Reynolds  number  (Res  =  4.1  x  104)  for  the  loca¬ 
tion,  extent  and  magnitude  of  the  pressure  plateau. 
Also,  the  shape  of  the  experimental  pressure  plateau 
shows  little  variation  for  Res  <  6.8  x  104,  thus  sug¬ 
gesting  that  the  computed  pressure  plateau  region  (for 
Res  =  2  x  104)  is  accurate.  The  computed  recovery  of 
the  surface  pressure  is  more  rapid  than  in  the  experi¬ 
ment,  however. 

The  computed  and  experimental  mean  skin  friction 
coefficient  Cf  —rwj \ PoqJJ ^  are  shown  in  Fig.  42.  The 
skin  friction  coefficient  increases  rapidly  downstream 
of  attachment.  The  computed  results  at  Res  =  2  x  104 
are  in  close  agreement  with  the  experimental  data  at 
Res  =  8.0  x  104  and  1.94  x  105  in  the  region  down¬ 
stream  of  reattachment. 


Fig.  41  Surface  pressure 


Fig.  42  Skin  friction  coefficient 

Zheltovodov  and  Schiilein13  has  correlated  experi¬ 
mental  data  for  the  scaled  separation  length  Lsep/Lc 
defined  by  (1)  for  the  expansion-compression  corner 
where  5e  is  the  incoming  boundary  layer  thickness 
(upstream  of  the  expansion  corner),  p2  is  the  pres¬ 
sure  after  the  shock  in  inviscid  flow,  ppi  is  the  plateau 
pressure  from  the  empirical  formula  of  Zukoski15 
Ppi  —Pe(\Me  +  1)  where  pe  and  Me  are  the  static 
pressure  and  freestream  Mach  number  upstream  of  the 
compression  corner  and  downstream  of  the  expansion 
fan.  In  the  computation,  the  location  is  taken  to  be 
x  =  65.  The  values  of  Me  and  p2  have  been  com¬ 
puted  using  inviscid  theory.  Also,  Rese  =  1.8  x  104  for 
les  (Rese  =  peUe5e/pe,  where  pe,  Ue  and  pe  are  com¬ 
puted  using  inviscid  theory).  The  experimental  data 
correlation  and  the  computed  result43  for  the  scaled 
separation  length  is  shown  in  Fig.  43.  The  computed 
value  is  consistent  with  a  linear  extrapolation  of  the 
experimental  data. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  les  accurately  reproduces  the  mean  surface 
pressure  and  rapid  recovery  of  the  skin  friction 
downstream  of  attachment  to  a  level  close  to  the 
undisturbed  upstream  in  agreement  with  experi¬ 
ment 


Fig.  43  Separation  length 

•  The  LES  predicts  a  mean  separation  length  con¬ 
sistent  with  a  linear  extrapolation  of  the  experi¬ 
mental  data 

2-D  Shock  Impingement 

The  shock  impingement  (Fig.  44)  is  analogous  to  the 
2-D  compression  corner  (Fig.  1).  An  incident  shock, 
generated  by  a  wedge  of  angle  a,  interacts  with  an 
equilibrium  supersonic  turbulent  boundary  layer.  For 
sufficiently  large  pressure  rise,  the  boundary  layer  sep¬ 
arates  and  a  system  of  shock  waves  and  expansions 
form.  The  mean  surface  pressure  distribution  is  simi¬ 
lar  to  the  2-D  compression  corner  at  the  same  overall 
pressure  rise.12  Reviews  are  presented  in  Green11  and 
Delery.12 


Fig.  44  2-D  shock  impingement 


Case  1:  Mach  2.3,  a  =  8° 

Gamier  et  a/50,51  simulated  a  Mach  2.3  impinging 
shock  at  Res  =  6.0  x  104.  The  incident  shock  angle 
is  32.4°  corresponding  to  a  wedge  angle  a  =  8°.  The 
wall  is  adiabatic.  Experimental  data  of  Deleuze52  and 
and  Laurent53  is  available.  The  flow  conditions  are 
summarized  in  Table  13.  The  overall  inviscid  pressure 
ratio  across  the  shocks  p2/pi  =  2.47. 

Table  13  Flow  Conditions 

Reference  Data  Res  x  10“ 4 

Gamier  et  al50,51  les  2.3  6.0 

Deleuze,52  Laurent53  E  2.3  6.0 

Details  of  the  computation  are  presented  in  Ta¬ 
ble  14.  A  structured  grid  LES  algorithm  was  uti¬ 
lized.  Three  computations  (denoted  by  LES-x  (a) 


to  LES-x  (c))  were  performed  using  the  Mixed  Scale 
Model,54-56  and  one  computation  using  miles  (de¬ 
noted  by  les-m  (d)).  The  inflow  conditions  were 
obtained  from  a  separate  computation  of  a  flat  plate 
adiabatic  boundary  layer.  The  inflow  velocity  satis¬ 
fies  the  compressible  Law  of  the  Wall  as  indicated  in 
Fig.  45  where  hwa  indicates  the  experimental  data 
obtained  using  Hot  Wire  Anemometry. 

Table  14  Details  of  Computations 


Reference 

Type 

Grid 

SA 

TA 

Gamier50,51 

LES-X  (A) 

255  x  55  x  151 

2nd 

2nd 

LES-X  (B) 

255  x  110  x  151 

2nd 

2nd 

LES-X  (C) 

510  x  55  x  151 

2nd 

2nd 

LES-M  (D) 

510  x  55  x  151 

2nd 

2nd 

Fig.  45  Velocity  (a  — ,  B  A,  c - ,  D  o,  hwa  o) 

The  mean  skin  friction  coefficient  is  shown  in  Fig.  46. 
The  computations  using  the  Mixed  Scale  Model  are 
in  close  agreement  with  experiment  and  each  other. 
The  miles  computation  (Case  D)  also  displays  close 
agreement  with  experiment,  and  exhibits  a  smaller 
separation  region  (consistent  with  higher  inflow  tur¬ 
bulence  level  near  the  wall50,51). 


250  300  350  400 

X 


Fig.  46  Skin  friction  (a  — ,  B  A,  c - ,  D  o,  hwa  o) 

The  evolution  of  the  compressible  displacement 
thickness  Si  is  displayed  in  Fig.  47.  All  computations 


shown  good  agreement  with  experiment  within  the 
interaction  region,  although  overestimating  in  the 
recovery  region  downstream  by  approximately  10%. 


Fig.  47  Displacement  thickness  (a  — ,  B  A,  c - , 

D  O,  HWA  o) 

The  mean  streamwise  velocity  profiles  at  x  = 
351  mm  (computations)  and  x  =  345  mm  (exper¬ 
iment)  are  shown  in  Fig.  48  where  lda  indicates 
the  experimental  measurements  using  Laser  Doppler 
Anemometer.  An  additional  profile  for  Case  A  is 
shown  at  x  =  345  mm.  The  x  location  is  imme¬ 
diately  downstream  of  reattachment  (Fig.  46).  Note 
the  velocity  is  the  ordinate  and  the  distance  from 
the  wall  (normalized  by  the  incoming  boundary  layer 
thickness)  is  the  abscissa.  The  offset  between  the  x  lo¬ 
cations  for  the  computations  and  experiment  accounts 
for  the  apparent  displacement  of  the  incident  shock 
in  the  computations.  The  overall  agreement  between 
computation  and  experiment  is  good. 


Fig.  48  Velocity  (a  ( x  —  351  mm)  — ,  a  (x  — 
345  mm)  —  •  — ,  B  A,  c - ,  d  o,  hwa  o,  lda  q) 

The  rms  streamwise  fluctuating  velocity  profiles  at 
x  =  326  mm  (computations)  and  x  =  320  mm  (ex¬ 
periment)  are  shown  in  Fig.  49,  and  at  x  =  351  mm 
(computations)  and  x  =  345  mm  (experiment)  in 
Fig.  50.  An  additional  profile  for  Case  A  is  shown 
at  x  =  320  mm  and  x  =  345  mm  in  Figs.  49  and 
50,  respectively.  The  profiles  in  Fig.  49  correspond 


to  the  separated  region,  and  the  profiles  in  Fig.  50  are 
located  downstream  of  the  reattachment  of  the  bound¬ 
ary  layer  (Fig.  46).  In  Fig.  49,  the  computed  profiles 
display  general  agreement  with  experiment  (lda).  In 
Fig.  50,  there  is  a  significant  difference  between  the 
hwa  and  lda  measurements  especially  in  the  middle 
of  the  boundary  layer.  The  computations  show  good 
agreement  with  the  lda  data. 


Fig.  49  Rms  velocity  (a  ( x  —  326  mm)  — ,  A  (x  — 
320  mm)  —  •  — ,  B  A,  c - ,  d  o,  lda  q) 


Fig.  50  Rms  velocity  (a  ( x  —  351  mm)  — ,  a  (x  — 

345  mm)  —  •  — ,  B  A,  c - ,  d  o,  hwa  o,  lda  q) 

In  summary,  the  following  conclusions  can  be  made: 

•  The  les  accurately  predicts  the  displacment 
thickness,  skin  friction  coefficient,  and  velocity 
profiles 

•  The  les  shows  reasonable  agreement  with  the 
LDA  measurements  of  rms  streamwise  turbulent 
fluctuations 

•  The  differences  between  the  Mixed  Model  and 
MILES  results  are  small 


3-D  Single  Fin 

The  3-D  single  fin  (Fig.  51)  is  a  standard  CFD 
validation  case  for  shock  wave  turbulent  boundary 
layer  interaction.  An  equilibrium  supersonic  turbulent 
boundary  layer  approaches  a  semi-infinite  fin  of  angle 
a.  The  deflection  of  the  flow  by  the  fin  generates  a 
swept  shock  system.  For  a  sufficiently  large  pressure 
rise,  the  boundary  layer  separates  and  a  A-shock  forms. 
The  flowfield  structure  is  well  understood.3,57-65 

Fin 


Fig.  51  3-D  single  fin 

Zheltovodov  et  a/57,58,60  identified  six  regimes  for 
the  3-D  single  fin  flowfield  depending  on  the  strength 
of  the  shock  wave.  Fig.  52  displays  the  surface  stream¬ 
line  patterns  for  each  of  the  regimeslx.  In  Regime  I, 
the  boundary  layer  is  unseparated  and  no  convergence 
of  surface  streamlines  is  observed.  In  Regime  II,  the 
streamlines  turn  approximately  parallel  to  the  invis- 
cid  shock  but  do  not  form  a  line  of  coalescence.  In 
Region  III,  a  primary  separation  line  (Si)  forms  corre¬ 
sponding  to  the  coalescence  of  the  surface  streamlines. 
An  attachment  line  (Ri)  forms  near  the  fin-plate  junc¬ 
ture.  A  line  of  secondary  separation  (52)  appears  lo¬ 
cated  between  the  primary  separation  (Si)  and  attach¬ 
ment  (Ri)  lines.  In  Regimes  IV  to  VI,  the  secondary 
separation  (52)  and  attachment  (R2)  lines  disappear 
and  finally  reappear. 

General  RANS  models  are  capable  of  predicting  the 
primary  features  of  the  flowfield  including  the  wall 
pressure  distribution  (albeit  with  some  notable  er¬ 
rors  for  strong  interactions),  and  location  of  separa¬ 
tion  and  attachment  lines.  However,  the  skin  friction 
distribution  is  not  accurately  predictedx  for  strong 
interactions.1,3,62  Moreover,  the  phenomenon  of 
the  appearance/disappearance/re- appearance  of  sec¬ 
ondary  separation  has  not  been  previously  predicted.61 

Cases  1-3:  Mach  3,  a  =  15°  and  Mach  4,  a  =  20° 
and  30.6° 

Thivet66  simulated  three  different  configurations  of 
the  3-D  single  fin.  The  flow  conditions  are  listed  in 
Table  15.  Experimental  data  of  Zheltovodov  et  a/57,58 

1XSW  identifies  the  trace  of  the  inviscid  shock. 

XA  notable  exception  is  the  prediction  of  the  modified 
Baldwin-Lomax  model  by  Panaras.2 


Fig.  52  3-D  single  fin  regimes 

and  Kim  et  al 67  is  available  for  Cases  1,3  and  2,  re¬ 
spectively.  The  computations  were  performed  at  the 
same  flow  conditions  as  the  experiment.  The  inviscid 
pressure  ratio  P2/P1  is  2.82,  5.21  and  9.5  for  Cases  1 
to  3,  respectively. 

Table  15  Flow  Conditions 


Case 

Reference 

Data 

Moo 

a 

Re§ 

xl0“ 

1 

Thivet66 

RANS 

3.0 

15° 

1.9 

Zheltovodov57’58 

E 

3.0 

15° 

1.9 

2 

Thivet66 

RANS 

4.0 

20° 

2.1 

Kim  et  al 67 

E 

4.0 

20° 

2.1 

3 

Thivet66 

RANS 

4.0 

30.6° 

1.6 

Zheltovodov  et  al 57 

E 

4.0 

30.6° 

1.6 

Details  of  the  computation  are  presented  in  Ta¬ 
ble  16.  Two  different  turbulence  models  were  exam¬ 
ined,  namely,  the  linear  Wilcox  k  —  oo  model68  (WI) 
and  a  weakly  nonlinear  Wilcox  k—co  model69  (WD+). 
The  WD+  model  is  an  extension  to  compressible  flows 
of  the  nonlinear  correction  of  Durbin.70  In  the  WD+ 
model,  the  coefficient  c ^  in  the  turbulent  eddy  viscos¬ 
ity  P't  —  is  defined  as  =  min(c°,  Jc^/s) 

where  s  =  S/uj ,  S  =  y/2SijSij  -  %SkkSkk  and  Sij  = 
1  (dUi / dxj  +  dUj  / dxi) .  Computations  were  performed 
using  the  GASPex  code.71  Inviscid  fluxes  are  com¬ 
puted  to  third  order  using  the  Roe  scheme  with 
MUSCL  reconstruction  and  a  min-mod  limiter.  Vis¬ 
cous  and  heat  transfer  terms  are  computed  using  sec¬ 
ond  order  central  differences.  A  typical  grid  is  shown 
in  Fig.  53. 

The  surface  streamlines  are  displayed  in  Figs.  54 
(WI),  55  (experiment)  and  56  (WD+)  for  Case  1,  and 
in  Figs.  57  (WI),  58  (experiment)  and  59  (WD+)  for 


Table  16  Details  of  Computations 


Case 

Reference 

Type 

Grid 

SA 

TA 

1 

Thivet66 

RANS 

128  x  160  x  160 

2nd 

ss 

2 

Thivet66 

RANS 

128  x  160  x  160 

2nd 

ss 

3 

Thivet66 

RANS 

112  x  160  x  160 

2nd 

ss 

Fig.  53  Grid  for  3-D  single  fin 


Case  3.  The  angle  of  the  primary  separation  line 
Si  is  underestimated  by  the  linear  (WI)  model  by 
1.4°  to  4°,  and  overestimated  by  the  weakly  nonlin¬ 
ear  (WD+)  model  by  1.4°  to  2.5°.  The  angle  of  the 
primary  attachment  line  R\  is  accurately  predicted 
by  both  models.  In  Case  1  (Regime  III),  the  weak 
secondary  separation  S2  observed  in  the  experiment 
(Fig.  55)  is  absent  in  in  the  computations  with  WI 
model  (Fig.  54).  For  the  WD+  model  (Fig.  56),  the 
distinct  changes  in  the  curvature  of  surface  stream¬ 
lines  occuring  midway  between  the  fin  and  primary 
separation  line  are  similar  to  incipient  secondary  sepa¬ 
ration  conditions  which  occurred  in  the  experiment  at 
a  lower  fin  angle.  In  Case  3  (Regime  VI),  the  strong 
secondary  separation  S2  and  attachment  R2  lines  ob¬ 
served  in  the  experiment  (Fig.  58)  are  evident  in  the 
WD+  model  (Fig.  59)  but  are  entirely  absent  in  the 
WI  model  (Fig.  57).  This  is  a  significant  achievement 
of  the  WD+  model. 

The  computed  and  experimental  pressure  distribu¬ 
tions  are  displayed  in  Figs.  60  to  62.  The  WD+  model 
shows  a  significant  improvement  compared  to  the  WI 
model,  and,  in  particular,  accurately  predicts  the  ap¬ 
pearance  of  the  secondary  shock  wave  associated  with 
the  secondary  separation. 

The  computed  and  experimental  skin  friction  co¬ 
efficient  for  Case  2  is  shown  in  Fig.  63.  The  peak 
skin  friction  for  the  WD+  model  is  reduced  by  35% 
compared  to  the  WI  model.  Two  series  of  measure¬ 
ments  are  shown,  where  the  solid  symbols  •  correspond 
to  corrected  values  of  the  skin  friction.72  The  WD+ 
model  displays  significantly  better  agreement  with  ex¬ 
periment  than  the  WI  model. 

The  improvement  achieved  by  the  WD+  model  is 


Fig.  54  Surface  streamlines  for  Case  1  (WI) 


Fig.  55  Surface  streamlines  for  Case  1  (exper) 


Fig.  56  Surface  streamlines  for  Case  1  (WD+) 


Cf  x  1000  (R  =  88.9  mm) 


Fig.  63  Skin  friction  for  Case  2 

associated  with  a  significant  reduction  in  the  peak  k 
within  the  primary  separation  vortex  core  and  in  the 
near  wall  flow  penetrating  to  the  place  of  formation  of 
the  secondary  separation  line.  Figs.  64  and  65  display 
computed  turbulence  kinetic  energy  contours  for  the 
WI  and  WD+  models,  respectively.  The  WD+  model 
reduces  the  peak  k  within  the  core  by  more  than  a 
factor  of  two,  thereby  increasing  the  size  of  the  inter¬ 
action.  These  results  also  support  the  conclusion  of 
Zheltovodov  et  a/57, 58,60  that  the  turbulence  amplifi¬ 
cation  in  the  reversed  flow  may  be  by  the  reason  for 
the  disappearance  of  the  secondary  separation  in  the 
regime  IV  -  V  (Fig.  52). 


Fig.  64  Turbulence  kinetic  energy  for  Case  3  (WI) 


Fig.  65  Turbulence  kinetic  energy  for  Case  3 
(WD+) 


In  summary,  the  following  conclusions  can  be  made: 

•  The  weakly  nonlinear  Wilcox-Durbin  (WD+) 
model  provides  a  significant  improvement  in  pre¬ 
diction  of  the  3-D  single  fin  interaction  in  compar¬ 
ison  with  the  Wilcox  (WI)  model.  In  particular, 

-  The  secondary  separation/attachment  and 
secondary  shock  are  captured  for  the  Mach  4, 
a  =  30.6°  configuration  (Regime  VI).  The 
linear  Wilcox  model  does  not  predict  sec¬ 
ondary  separation/attachment  and  the  sec¬ 
ondary  shock  for  this  case. 

—  The  surface  pressure  distribution  is  in  close 
agreement  with  experiment  for  all  three 
cases,  except  for  the  position  of  upstream  in¬ 
fluence  (i.e.,  the  location  of  initial  pressure 
rise)  for  the  strongest  interaction 

•  The  improvement  in  prediction  by  the  Wilcox- 
Durbin  model  is  associated  with  a  reduction  in 
the  computed  turbulence  kinetic  energy  k  in  the 
vortex  core 

3-D  Double  Fin 

The  3-D  double  fin  (crossing  shock  wave)  (Fig.  66) 
has  become  a  canonical  validation  case  for  CFD  mod¬ 
eling  of  shock  wave  turbulent  boundary  layer  interac¬ 
tions.  An  equilibrium  supersonic  turbulent  boundary 
layer  approaches  a  channel  defined  by  two  semi-infinite 
fins  of  angles  a±  and  a2  separated  by  a  distance  We 
at  the  entrance.  Additional  geometric  parameters  in¬ 
clude  the  minimum  (throat)  width  Wt  and  the  offset 
(if  any)  of  the  midpoint  of  the  throat  relative  to  the 
midpoint  of  the  entrance.  The  flowfield  structure  of 
the  double  fin  interaction  has  been  investigated  by  sev¬ 
eral  researchers73-91  and  the  basic  topological  features 
have  been  accurately  predicted  for  weak  and  moder¬ 
ate  interactions;  however,  a  complete  understanding 
of  the  flow  has  not  yet  been  achieved.  An  example  of 
the  mean  streamline  structure  is  shown  in  Fig.  67  from 
Thivet  et  al 69  for  a  Mach  4,  7°  x  11°  double  fin  inter¬ 
action.  The  principal  features  include  the  separated 
boundary  layer,  vortex,  and  entrained  flow.  The  com¬ 
puted  surface  heat  transfer  coefficient  contours  and 
surface  streamlines  are  shown  in  Fig.  68  and  the  cor¬ 
responding  experimental  surface  streamline  pattern  in 
Fig.  69.  The  primary  separation  lines  (denoted  S± 
and  S2  in  Fig.  69)  are  evident  in  both  computation 
and  experiment,  and  correspond  to  the  liftoff  of  the 
streamlines  from  the  surface  in  Fig.  67.  The  attach¬ 
ment  lines  (denoted  R\  and  R2  in  Fig.  69)  correspond 
to  the  impingement  of  the  streamlines  which  demar¬ 
cate  the  boundaries  of  the  counter-rotating  vortex  pair 
formed  by  the  separated  flow.  Secondary  separation 
lines  (53  and  S4  in  Fig.  69)  and  attachment  line  (de¬ 
noted  R3  in  Fig.  69)  are  also  evident. 


Fin 


Fig.  66  3-D  double  fin 


Fig.  67  3-D  double  fin  flow  structure 


Ch  x  1 03  :  0.67  1.00  2.00  2.70  3.00  4.00 


Fig.  68  Computed  surface  heat  transfer 


Fig.  69  Surface  flow  visualization 

Cases  1-2:  Mach  4,  a\  —  7°, 02  =  11°  and  Mach  4, 

on  =  15°,  a2  =  15° 

Thivet69  simulated  two  configurations  of  the  3-D 
double  finxl.  Experimental  data  of  Zheltovodov  et 

X1A  third  configuration  (Moo  =  4,  a\  —  7°,  <22  =  7°)  is  not 
presented  here,  as  it  is  a  weaker  interaction.  It  is  accurately 


al77~80  is  available  for  both  cases.  The  flow  condi¬ 
tions  are  listed  in  Table  17.  The  computations  were 
performed  at  the  same  flow  conditions  as  the  exper¬ 
iment.  The  experimental  configuration  for  Case  1  is 
shown  in  Fig.  70. 

Table  17  Flow  Conditions 


Case 

Reference 

Data 

Moo 

Oil 

OL2 

Res 

xl0“ 

1 

Thivet69 

RANS 

3.92 

70 

11° 

3.0 

Zheltovodov  et  al77~80 

E 

3.92 

7° 

11° 

3.0 

2 

Thivet69 

RANS 

3.92 

15° 

15° 

3.0 

Zheltovodov  et  al77~80 

E 

3.92 

15° 

15° 

3.0 

Fig.  70  3-D  double  fin  configuration 

Details  of  the  computations'*11  are  presented  in  Ta¬ 
ble  18.  Five  different  turbulence  models  were  consid¬ 
ered,  namely,  the  (linear)  Wilcox  k—oo  (WI)  model,68 
and  four  variants  of  the  Wilcox  model  based  on  the 
requirement  for  weak  realizability92  of  the  computed 
Reynolds  stressesxm.  The  variants  are  the  Wilcox- 
Moore  (WM),  Wilcox-Durbin  (WD)  and  two  modi¬ 
fications  thereof  (WM+  and  WD+).  Details  of  the 
models  are  presented  in  Thivet  et  al.m  A  detailed 
grid  refinement  study69  was  performed  to  quantify  the 
uncertainty  in  the  computatations. 

Table  18  Details  of  Computations 


Case 

Reference 

Type 

Grid 

SA 

TA 

1 

Thivet69 

RANS 

96  x  80  x  160 

2nd 

SS 

2 

Thivet69 

RANS 

192  x  200  x  104 

2nd 

ss 

The  computed  and  experimental  surface  pressure 
and  heat  transfer  coefficientxlv  on  the  Throat  Middle 
Line  (TML)XV  for  Case  1  are  shown  in  Fig.  71.  The 

predicted  by  general  rans  models. 
xllData  are  shown  for  the  finest  grid. 

xm Realizability  implies  that  the  normal  kinematic  Reynolds 
stresses  v!u\  v'v'  and  w'w'  are  non-negative,  and  the  Reynolds 
shear  stresses  u'v'  ,u'w'  and  v'w'  are  bounded  by  the  Schwartz 
inequality,  e.g .,  (u'v')2  <  u'u'  v'v' . 

xlvThe  heat  transfer  coefficient  =  qw/ Poo UooCp(Tw  —  Taw) 
where  Taw  is  the  adiabatic  wall  temperature. 

xvThe  Throat  Middle  Line  is  the  streamwise  datum  through 
the  center  of  the  throat. 


surface  pressure  is  accurately  predicted  by  all  models 
with  only  slight  differences  in  the  plateau  pressure  at 
x  =  80  —  100  mm.  Similar  good  agreement  with  exper¬ 
iment  and  insensitivity  of  the  predictions  to  the  tur¬ 
bulence  model  are  observed  in  the  spanwise  pressure 
distributions  at  three  streamwise  locations  (Fig.  72). 
This  is  consistent  with  previous  results1  which  showed 
a  marked  insensitivity  of  the  computed  surface  pres¬ 
sure  to  the  turbulence  model  employed. 

The  computed  surface  heat  transfer  coefficient  Ch  on 
the  TML  (Fig.  71)  and  at  three  streamwise  locations 
(Fig.  73)  significantly  overpredicts  the  peak  experi¬ 
mental  value  for  all  models.  This  is  similar  to  previous 
results  using  other  RANS  models  which  did  not  incor¬ 
porate  a  specific  realizability  effect.1 


Fig.  71  Surface  pressure  and  heat  transfer  (ai  = 

7°,  a2  =  11°) 


Fig.  72  Surface  pressure  (ai  =  7°,a2  =  11°) 

The  computed  and  experimental  surface  pressure 
and  heat  transfer  coefficient  on  the  Throat  Middle  Line 
(TML)  for  Case  2  are  shown  in  Fig.  74  for  the  WI  and 
WM+  models.  All  three  models  underestimate  the  lo¬ 
cation  of  the  pressure  rise  immediately  following  the 
plateau  and  overestimate  the  peak  pressure.  The  sub¬ 
sequent  expansion  (associated  with  the  corner  of  the 


Fig.  73  Surface  pressure  (ai  =  7°,a2  =  11°) 

fins)  is  more  accurately  predicted  by  the  WI  model. 
The  computed  peak  heat  transfer  significantly  overes¬ 
timates  the  experiment  for  all  models. 


Fig.  74  Surface  pressure  and  heat  transfer  (ai  = 

15°,  a2  =  15°) 

Thivet  examined  the  effect  of  imposing  the  turbu¬ 
lence  length  scale  limitation  proposed  by  Coakley  and 
Huang.93  The  turbulent  eddy  viscosity  fit  =  pcpk^^l 
and  the  turbulence  length  scale  l  is  defined  by 

l  =  min  [l^k' '*!«,] 


where 


^°g  -  3/4 


ell  ~d 


where  k  =  0.41  is  von  Karman’s  constant  and  d 
is  the  distance  normal  to  the  wall.  The  dissipation 
term  in  the  transport  equation  for  k  is  modeled  using 
c o  =  k}/2t.  Fig.  75  displays  the  ratio  of  the  computed 
turbulence  length  scale  l  (without  the  Coakley  and 
Huang  correction)  to  the  logarithmic  expression  £\og 
for  Case  1  at  x  =  46  mm  for  six  spanwise  locations. 
The  transverse  location  zr  is  zero  at  the  fin  and  unity 


at  the  location  corresponding  to  the  inviscid  shock  gen¬ 
erated  by  the  11°  fin.  It  is  evident  that  the  ratio  does 
not  significantly  exceed  unity,  thereby  implying  that 
the  Coakley  and  Huang  correction  would  have  little 
effect  on  the  computed  solution. 

N=  0  1  23  45 


l\g(+N) 

Fig.  75  Turbulence  length  scale  ratio  (ai  =  7°,  02  = 
11°) 

Thivet  also  examined  the  effect  of  incorporating  a 
compressibility  correction68  in  the  k— 00  modelxvl.  The 
threshold  value  of  the  turbulence  Mach  number  Mt  = 
^jTkfa  (where  a  is  the  local  speed  of  sound)  is  0.25 
for  the  compressibility  correction.  Below  this  value, 
the  correction  is  inoperative.  The  maximum  value  of 
Mi  using  the  WI  model  is  0.3  yielding  an  increase  of 
only  4%  in  the  dissipation  term,  thereby  implying  that 
this  correction  (if  implemented)  would  have  negligible 
effect  on  the  computation  using  the  WI  model.  The 
values  of  Mt  for  the  WD,  WD+  and  WM+  models 
are  everywhere  below  the  threshold,  and  therefore  the 
compressibility  correction  (if  implemented)  would  have 
no  effect  on  these  results. 

In  summary,  the  following  conclusions  can  be  made: 

•  The  computed  surface  pressure  is  accurately  pre¬ 
dicted  by  the  five  different  turbulence  models  for 
Case  1  (ai  =  7°,  a2  =  11°).  Some  discrepancies 
between  computed  and  experimental  surface  pres¬ 
sure  are  evident  for  Case  2  (a±  =  15°,  012  =  15°). 

•  The  peak  heat  transfer  is  not  accurately  predicted 
by  any  of  the  models  for  either  case. 

•  The  proposed  length  scale  correction  of  Coakley 
and  Huang,93  and  the  compressibility  modifica¬ 
tions  of  Zeman94  or  Sarkar,95  would  not  signifi¬ 
cantly  change  the  computed  flowfields. 

xvlThe  compressibility  correction  concept,  introduced  by  Ze¬ 
man94  and  Sarkar95  to  account  for  the  observed  decrease  in 
the  spreading  rate  of  compressible  mixing  layers,  leads  to  an 
undesirable  reduction  in  skin  friction  for  turbulent  boundary 
layers.68,93,96 


Case  3:  Mach  5,  ai  =  02  =  18° 

Schmisseur  and  Gaitonde97  and  Panaras98  simu¬ 
lated  a  Mach  5,  a±  =  «2  =  18°  double  fin  at  Res  = 
1.4  x  105  and  isothermal  wall  conditions  ( Tw/Taw  = 
0.76).  Experimental  data  of  Schuelein  and  Zhel- 
tovodov99,100  (an  additional  description  is  presented 
in  Zheltovodov  and  Maksimov101)  is  available  for  this 
case.  The  flow  conditions  are  shown  in  Table  19. 
The  computations  were  performed  at  the  same  flow 
conditions  as  the  experiment.  The  experimental  con¬ 
figuration  is  shown  in  Fig.  76. 

Table  19  Flow  Conditions 


Reference 

Data 

Moo 

ai 

Q!2 

Res 

Schmisseur  et  al 97 

RANS 

4.96 

18° 

18° 

1.39  x  10J 

Panaras98 

RANS 

4.96 

18° 

18° 

1.39  x  10J 

Schiilein  and 

E 

4.96 

18° 

18° 

1.39  x  10J 

Zheltovodov99  101 


Fig.  76  Double  fin  configuration 

Details  of  the  computations  are  presented  in  Ta¬ 
ble  20.  Schmisseur  and  Gaitonde  utilized  the  k  —  e 
model102,103  with  low  Reynolds  number  and  com¬ 
pressibility  corrections.  Panaras  employed  a  modified 
Baldwin-Lomax  model.2 

Table  20  Details  of  Computations 


Reference 

Type 

Grid 

SA 

TA 

Schmisseur  et  al97 

RANS 

153  x  143  x  103 

2nd 

SS 

Panaras98 

RANS 

127  x  101  x  91 

2nd 

ss 

The  experimental  surface  oil  flow  is  presented  in 
Fig.  77,  and  the  computed  surface  streamlines  in 
Fig.  78  (Schmisseur  and  Gaitonde)  and  Fig.  79  (Pa¬ 
naras).  The  experiment  displays  a  complex  structure 
of  separation  (5)  and  attachment  (R)  lines.  The  initial 
separation  lines  (Si  and  52)  generated  by  the  indi¬ 
vidual  fin  shock  boundary  layer  interactions  coalesce 
to  form  a  node  at  the  centerline.  Corresponding  at¬ 
tachment  lines  (R±  and  R2)  are  located  near  the  fins. 
Additional  separation  (53  to  S$)  and  attachment  lines 
(Rs  to  Rq)  form  downstream.  The  central  point  (51) 
is  a  saddle  point  in  the  experiment.  In  the  same  re¬ 
gion,  the  computation  of  Schmisseur  and  Gaitonde 
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(Fig.  78)  displays  a  node  (iVi),  while  the  computa¬ 
tion  of  Panaras  (Fig.  79)  appears  in  closer  agreement 
with  experiment.  Nonetheless,  both  computed  surface 
streamline  plots  show  overall  good  agreement  with  ex¬ 
periment. 


Fig.  77  Experimental  surface  oil  visualization 

spanwise  pressure  survey  locations 


Fig.  78  Surface  streamlines 


Fig.  79  Surface  streamlines 

The  computed  and  experimental  surface  pressure 
on  the  centerline  is  displayed  in  Fig.  80  (Schmisseur 
and  Gaitonde)  and  Fig.  81  (Panaras).  The  experimen¬ 
tal  pressure  profile  displays  a  rise  due  to  the  crossing 
shock  interaction,  and  a  subsequent  drop  beginning  at 
x  =  42 5  due  to  the  expansion  originating  from  the 
corners  C  in  Fig.  76.  The  subsequent  pressure  rise 
is  associated  with  the  reflection  of  the  shocks  from 
the  fin  surfaces.  The  computed  surface  pressure  of 


Schmisseur  and  Gaitonde  shows  significant  deviation 
from  experiment  beginning  at  the  initial  plateau  in 
surface  pressure.  The  computed  surface  pressure  of 
Panaras  shows  good  agreement  up  to  the  beginning 
of  the  decrease  in  pressure  at  x/S  =  42.  However,  a 
significant  discrepancy  appears  thereafter. 

The  computed  and  experimental  surface  pressure  at 
three  streamwise  locations  are  shown  in  Figs.  82  to  87. 
The  results  of  Schmisseur  and  Gaitonde  are  presented 
in  Figs.  82,  84  and  86,  and  the  results  of  Panaras  in 
Figs.  83,  85  and  87.  Both  computations  show  good 
agreement  with  experiment  at  x  =  26.65  (Figs.  82 
and  83).  This  location  is  upstream  of  the  interaction 
of  the  crossing  shocks,  and  thus  represents  individual 
single  fin  interactions.  At  x  =  37.15  (located  at  ap¬ 
proximately  the  second  focal  point  N 2  as  indicated  in 
Fig.  78)  both  computed  profiles  show  good  agreement 
with  experiment,  except  for  a  overshoot  in  centerline 
pressure  in  Fig.  84  (Schmisseur  and  Gaitonde).  At 
x  =  52.95  (see  Fig.  78),  both  computations  show  gen¬ 
erally  good  agreement  with  experiment  except  for  an 
overshoot  in  pressure  at  the  centerline. 


Fig.  80  Surface  pressure  on  centerline 


Fig.  81  Surface  pressure  on  centerline 
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Fig.  82  Surface  pressure  at  x  —  26.66 


Fig.  84  Surface  pressure  at  x  —  37.1 6 


Fig.  83  Surface  pressure  at  x  =  26.6$ 

In  summary,  the  following  conclusions  can  be  made: 

•  The  computed  surface  streamlines  are  in  good 
agreement  with  experiment.  The  detailed  struc¬ 
ture  of  separation,  attachment  and  foci  are  gen¬ 
erally  accurately  reproduced  in  the  simulations. 

•  The  computed  surface  pressure  shows  general 
agreement  with  experiment  during  the  initial  por¬ 
tion  of  the  interaction  and  off  centerline  through¬ 
out  the  entire  interaction.  However,  significant 
differences  between  computed  and  experimental 
surface  pressure  are  evident  on  the  centerline  of 
the  interaction. 

Case  4:  Mach  5,  ai  —  a2  —  23° 

Schmisseur  and  Gaitonde82, 86,88,97  and  Panaras98 
simulated  a  Mach  5,  a\  =  «2  =  23°  double  fin 
at  Res  =  1.4  x  105  and  isothermal  wall  conditions 
(: Tw/Taw  =  0.76).  Experimental  data  of  Schulein  and 
Zheltovodov99,100  (an  additional  description  is  pre¬ 
sented  in  Zheltovodov  and  Maksimov101)  is  available 
for  this  case.  The  flow  conditions  are  shown  in  Ta¬ 
ble  21.  The  computations  were  performed  at  the  same 


Fig.  85  Surface  pressure  at  x  —  37.1 6 

flow  conditions  as  the  experiment.  The  experimental 
configuration  is  shown  in  Fig.  88. 

Table  21  Flow  Conditions 


Reference 

Data 

Moo 

Oil 

OL2 

Res 

Schmisseur  et  al97 

RANS 

4.96 

23° 

23° 

1.39  x  10J 

Panaras98 

RANS 

4.96 

23° 

23° 

1.39  x  10J 

Schulein  and 

E 

4.96 

23° 

23° 

1.39  x  10J 

Zheltovodov"  101 

Details  of  the  computations™1  are  presented  in  Ta¬ 
ble  22.  Schmisseur  and  Gaitonde  utilized  the  k  —  e 
model102,103  with  low  Reynolds  number  and  com¬ 
pressibility  corrections.  Panaras  employed  a  modified 
Baldwin-Lomax  model.2 

The  experimental  surface  oil  flow  is  presented  in 
Fig.  89,  and  the  computed  surface  streamlines  in 
Figs.  90  (Schmisseur  and  Gaitonde)  and  91  (Panaras). 
Good  agreement  is  observed  between  the  computations 

xvl1  Schmisseur  and  Gaitonde  performed  a  grid  refinement  study 
using  meshes  of  91  x  72  x  52,  181  x  143  x  103  and  229  x  179  x  129 
and  confirmed  that  the  finest  grid  solution  was  grid  converged. 


Table  22  Details  of  Computations 


Fig.  89  Experimental  surface  oil  visualization 


Fig.  87  Surface  pressure  at  x  —  52.9$ 


Fig.  88  Double  fin  configuration 


and  experiment. 


The  computed  and  experimental  surface  pressure  on 
the  centerline  is  displayed  in  Fig.  92  (Schmisseur  and 
Gaitonde)  and  Fig.  93  (Panaras).  The  experimental 
profile  is  similar  to  the  18°  x  18°  configuration  except 
for  the  absence  of  the  second  shock  reflection  on  the 
centerline  due  to  the  limited  length  of  the  flat  plate. 
The  computed  profile  by  Schmisseur  and  Gaitonde 
shows  good  agreement  with  experiment  up  to  the  peak 
experiment  pressure,  but  significantly  overestimates 
experiment  downstream.  The  computed  profile  by  Pa¬ 
naras  also  shows  good  agreement  with  experiment  up 
to  the  peak  pressure,  but  does  not  display  the  expan¬ 
sion  regionxvm  observed  in  the  experimentq. 

The  computed  and  experimental  surface  pressure 

XV111A  separate  coarse  grid  computation  was  performed  by  Pa¬ 
naras104  wherein  the  effect  of  the  finite  plate  length  was  ap¬ 
proximated  by  adding  a  divergent  section  to  the  flat  plate. 
The  computed  centerline  pressure  displayed  significantly  better 
agreement  with  experiment. 


Fig.  91  Surface  streamlines 


Fig.  92  Surface  pressure  on  centerline 


Fig.  93  Surface  pressure  on  centerline 

at  three  streamwise  locations  are  shown  in  Figs.  94 
to  99.  The  results  of  Schmisseur  and  Gaitonde  are 
shown  in  Figs.  94,  96  and  98,  and  the  results  of  Pa- 
naras  in  Figs.  95,  97  and  99.  Both  computations  are  in 
close  agreement  with  the  available  experimental  data 
at  x  =  22.15  which  is  upstream  of  the  crossing  shock 
interaction  (see  Figs.  92  and  93).  At  x  =  32.65,  lo¬ 


cated  in  the  midst  of  the  crossing  shock  interaction, 
both  computations  show  close  agreement  with  exper¬ 
iment,  although  somewhat  overestimating  the  center- 
line  pressure.  At  x  =  48.45,  the  agreement  between 
computation  and  experiment  is  good  except  in  the  im¬ 
mediate  vicinity  of  the  centerline  where  the  pressure 
is  overestimated. 


Fig.  95  Surface  pressure  at  x  —  22.15 

In  summary,  the  following  conclusions  can  be  made: 

•  The  computed  surface  streamlines  display  good 
agreement  with  experiment 

•  The  computed  surface  pressure  shows  general 
agreement  with  experiment  during  the  initial  por¬ 
tion  of  the  interaction  and  off  centerline  through¬ 
out  the  entire  interaction.  However,  significant 
differences  between  computed  and  experimental 
surface  pressure  are  evident  on  the  centerline,  sim¬ 
ilar  to  Case  3. 

Conclusions 

An  evaluation  of  the  CFD  capability  for  predic¬ 
tion  of  shock  wave  turbulent  boundary  layer  inter¬ 
action  was  performed  under  the  auspices  of  NATO 
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Fig.  96  Surface  pressure  at  x  —  32.66 


Fig.  97  Surface  pressure  at  x  —  32.6$ 

RTO  Working  Group  10.  Five  separate  configurations 
were  considered,  namely,  the  2-D  compression  corner, 
2-D  expansion-compression  corner,  2-D  shock  impinge¬ 
ment,  3-D  single  fin  and  3-D  double  fin.  Recent  DNS 
s,  LES  and  RANS  simulations  were  solicited.  The  com¬ 
puted  results  were  compared  with  experimental  data 
where  available.  The  following  conclusions  may  be 
drawn: 

•  DNS  and  LES 

-  DNS  and  LES  of  shock  wave  turbulent  bound¬ 
ary  layer  interaction  is  an  active  area  of 
research.  In  the  previous  review  of  shock 
wave  turbulent  boundary  layer  interaction 
by  Knight  and  Degrez1  published  in  1998, 
there  were  no  DNS  or  LES  results  available. 
The  present  paper  includes  several  DNS  and 
LES  results  for  2-D  shock  wave  turbulent 
boundary  layer  interactions. 

—  A  comprehensive  assessment  of  the  capabil¬ 
ity  of  DNS  and  LES  was  not  possible  for  two 
reasons.  First,  the  DNS  and  LES  simulations 


20 


Fig.  99  Surface  pressure  at  x  —  48.4$ 

have  been  performed  at  lower  Reynolds  num¬ 
bers  Res  than  the  experiment  for  all  cases  ex¬ 
cept  one  (i.e.y  the  2-D  shock  impingement). 
Second,  the  DNS  and  LES  have  been  per¬ 
formed  for  nominally  2-D  flows.  Both  rea¬ 
sons  are  attributable  to  the  computational 
cost  of  DNS  and  LES. 

-  Although  the  set  of  comparisons  between 
DNS  and  LES  with  experiment  was  lim¬ 
ited,  nevertheless  it  is  clear  that  significant 
progress  has  been  achieved  in  the  predic¬ 
tion  of  shock  wave  turbulent  boundary  layer 
interaction  using  DNS  and  LES.  In  several 
cases,  the  computations  were  performed  at 
Reynolds  numbers  close  to  the  experiment 
(and,  in  one  case,  equal  to  experiment),  and 
good  agreement  with  experimental  data  was 
achieved.  Detailed  comparisons  included 
surface  pressure,  skin  friction,  velocity,  tem¬ 
perature  and  turbulence  profiles.  However, 
some  significant  discrepancies  are  evident, 
due  possibly  to  the  differences  in  Reynolds 


number. 

—  No  dns  or  les  results  for  surface  heat  trans¬ 
fer  in  shock  wave  turbulent  boundary  layer 
interaction  were  available.  This  is  partic¬ 
ularly  disappointing,  since  prior  RANS  sim¬ 
ulations  of  strongly  separated  shock  wave 
turbulent  boundary  layer  interactions  have 
failed  to  accurately  predict  heat  transfer.1 

•  RANS 

-  New  ideas  in  modeling  have  reinvigorated 
RANS  simulations  for  shock  wave  turbulent 
boundary  layer  interaction.  These  new  ideas 
include  realizability  (weak  nonlinearity)  and 
specific  physical  models  which  incorporate 
known  flowfield  behavior. 

—  The  weakly  nonlinear  Wilcox-Durbin  model 
successfully  predicts  secondary  separation 
for  the  3-D  single  fin  in  agreement  with  ex¬ 
periment. 

—  The  linear  and  weakly  nonlinear  Wilcox- 
based  models  fail  to  accurately  predict  sur¬ 
face  heat  transfer  for  the  separated  3-D  dou¬ 
ble  fin. 

-  The  surface  streamline  pattern  for  the  3-D 
double  fin  Case  Nos.  3  and  4  is  predicted 
by  both  the  conventional  k  —  t  and  Baldwin- 
Lomax-Panaras  models. 

-  The  surface  pressure  for  the  3-D  double  fin  is 
accurately  predicted  by  the  conventional  k—t 
and  Baldwin-Lomax-Panaras  models  only  in 
the  initial  region  of  the  interaction.  Signif¬ 
icant  differences  between  the  computational 
and  experimental  centerline  pressure  are  ev¬ 
ident  downstream. 

Future  Needs 

The  following  research  needs  for  shock  wave  turbu¬ 
lent  boundary  layer  interaction  may  be  identified: 

•  dns  and  les  simulations  are  needed  at  Reynolds 
numbers  corresponding  to  the  experiment.  This  is 
essential  to  end  speculation  regarding  the  reasons 
for  discrepancy  between  simulation  and  experi¬ 
ment. 

•  dns  and  les  simulations  are  needed  for  3-D  con¬ 
figurations.  All  dns  and  les  results  to  date  have 
been  for  nominally  2-D  flows. 

•  dns  and  les  simulations  are  needed  to  compare 
with  experimental  data  for  heat  transfer.  Neither 
dns,  les  nor  rans  models  to  date  have  success¬ 
fully  predicted  surface  heat  transfer  in  strongly 
separated  shock  wave  turbulent  boundary  layer 
interactions. 


•  RANS  modeling  needs  continued  emphasis  on 
weakly  nonlinear  corrections  to  two  equa¬ 
tion  models,  and  continued  development  of 
configuration-specific  modifications  to  simple 
models. 
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